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Abstract 

The electronic properties of thin films present quantum-size effects, which 
are a consequence of the finite size of the system. Here we focus on the investi- 
gation of these effects on the electronic energy loss of charged particles moving 
parallel with thin metallic films. The energy loss is calculated, within linear- 
response theory, from the knowledge of the density-response function of the 
inhomogeneous electron system, which we evaluate either in the random-phase 
approximation or with the use of an adiabatic and local exchange-correlation 
kernel. 
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1 Introduction 

The interaction of moving ions with solids has represented an active field of basic and 
applied physics Jl], 0. Charged particles moving near metallic surfaces loose energy 
as a consequence of the creation of different kind of excitations in the metal, such as 
electron- hole pairs and both bulk and surface plasmons|3|, |J. Also, the theoretical 
understanding of the electronic excitations is well known to be relevant in surface 
physics, as these modes are invariably involved in a variety of surface spectroscopies [5- 
9]- 

Quantum-size effects (QSE) were first investigated by Schulte[]TD| within a jellium 



model of the thin film and later by Feibelman|TT|, showing that the various physical 
properties exhibit an oscillatory behaviour as a function of the thickness of thin metallic 
films. These effects, which decrease as the size of the thin film increases, can be 
observed experimentally . QSE on the surface energy and the work function have 



been examined recently, within a stabilized jellium model of the electron system [|13| . 

In this paper we focus on the investigation of QSE on the energy-loss spectra of 
charged particles moving parallel with metallic slabs. The energy loss is calculated, 
within linear response theory, from the knowledge of the density-response function of 
the inhomogeneous electron system, which we evaluate either in the random-phase 
approximation or with the use of an adiabatic and local exchange-correlation kernel. 
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In Section 2 we present general expressions for the energy loss of charged particles 
moving along a definite trajectory. The results of our self-consistent calculations are 
presented in Section 3, and in Section 4 our main conclusions are summarized. Atomic 
units are used throughout, i.e., e 2 = H = m e = 1. 



2 Theory 

We consider a recoiless particle of charge Z\ moving with constant velocity v along a 
definite trajectory at a fixed distance z from the planar surface of a bounded three- 
dimensional electron gas that is translationally invariant in two directions. The energy 
that the moving particle looses per unit path length, i.e., the so-called stopping power 
of the electron system may be obtained from the gradient along the particle trajectory 
of the self-consistent potential set up by the particle in the electron system and eval- 



uated at the position of the particle W\. Within linear response theory, one finds fT5 



dE 2 f dqii f°° 

-^ = - v zl ij2^L ^"^(*.*;%")*("-<id-v), (1) 

where qy is the momentum transfer in the plane of the surface, u represents the energy 
transfer, and W(z, z\ qy, lu) is the screened interaction 

W(z,z';q h u) = v(z,z';q\\) + J dzi J dz 2 v(z,z 1 ;q\\)x(zi,z 2 ;q\\,uj)v(z2,z';q\\), (2) 

v(z,z';q\\) and x{ z i z ''i <Z||> u ) being two-dimensional Fourier transforms of the bare 
Coulomb potential and the density-response function, respectively fl6[ . 



The stopping power of the electron system can be described by means of -P(w), the 
total probability for the moving particle to exchange energy uj with the medium: 

dE 1 r°° , _,, x /n . 

dwujP(uj), (3) 



dx v Jo 
where 



P(u) = ± dq x lmW(z,z;q h uj), (4) 
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with q\\ = \Jq% + {uj/v) 2 and q x being the momentum transfer along the particle tra- 
jectory. 

The key ingredient of our energy-loss calculations is the density-response function 
of the electron system, which is known to satisfy the following integral equation W71 

x(z,z';Q\\,u) = x°(z,z';q\\,u) + J dz x J dz 2 x°(z,z';q {l ,uj) 

x [v(z 1 ,z 2 ;q\\) + f ms {z 1 ,Z2]q\\,u)\ x( z 2, z'; q\\, u), (5) 

where x°{ z i z x \ Qh u ) ls the density-response function of non-interacting electrons mov- 
ing in the effective Kohn-Sham potential of density-functional theory (DFT)|18[, and 
the kernel f xc (z, z'\ q\\,uj) accounts for exchange-correlation (xc) effects beyond a time- 
dependent Hartree approximation. 



Exchange-correlation effects are usually introduced within the local-density ap- 
proximation (LDA) of DFT, by replacing the xc potential at z by that of a uniform 
electron gas of density n(z). The xc kernel entering Eq. (El) is then set either equal 
to zero [this is the random-phase approximation (RPA)] or equal to the static (u = 0) 
xc kernel 

f£ DA {z,z';q h u) 



dv xc (n) 



dn 



S(z-z'). (6) 

n=n(z) 



This is the so-called adiabatic local-density approximation (ALDA). 

We consider a jellium slab of thickness a normal to the z axis, consisting of a fixed 
uniform positive background of density: 

!n, —a < 
(7) 
0, elsewhere, 

plus a neutralizing cloud of interacting electrons of density n(z). The positive-background 
charge density is n = qp/3ir 2 , where qF = (97r/4) 1 / 3 /r< j is the Fermi momentum and 
r s is the Wigner-Seitz radius. 



To compute x{ z i z ''i <7||> u )-> we follow the method described in Ref. [[19]. We first 
assume that n(z) vanishes at a distance Zq from either jellium edge, and expand the 
one-electron wave functions in a Fourier sine series. The distance zq and the number 
of sine functions kept in the expansion of the wave functions are chosen sufficiently 
large for our calculations to be insensitive to the precise values employed. We then 
introduce a double-cosine representation for the density-response function, and find 
explicit expressions for the screened interaction and the energy-loss probability in 



terms of the Fourier coefficients of the density-response function [|T5 . 

QSE are originated in the quantization of the energy levels normal to the surface. 
As the slab-thickness a increases new subbands for the z motion become occupied, 
thereby leading to oscillatory functions of a with and amplitude that decays approx- 
imately linearly with a and a period that equals Xp = 2k /qp being the Fermi 
wavelength. These size effects are responsible for the oscillatory behaviour of the 
electronic density induced in the electron system by the external probe and also for 
the oscillations of the energy-loss function and the stopping power with the system 
size. The results presented below correspond to slabs with the number n of occupied 
subbands in the range n = 11 — 14, for which a ~ 4 — 7 \p- 



3 Results 

We have investigated thin-film effects on the energy-loss spectra of aluminum slabs, 
for which r s = 2.07 and = 6.77 a (a is the Bohr radius, ao = 0.529 A). We set 
Z\ — ±1 and our results can then be used for arbitrary values of Z\, as the energy- loss 
probability is, within linear- response theory, proportional to Z\. 

The main ingredient of our calculations is the energy-loss function ImW^;?, z\ q\\,ou). 
Fig. 1 shows this quantity, as obtained from Eq. (0) with use of either the RPA or 
the ALDA density-response function, versus the slab-thickness a and for z = 2a , 
q\\ = OAqp, and uo = 0.2 lo p [uo p = (Ann) 1 / 2 is the classical plasma frequency of a 
uniform electron gas of density n]. As a increases, new subbands for the z motion 



become occupied. For r s = 2.07, the n = 11, 12, 13, 14 subbands fall below the Fermi 
level for a n = 4.95, 5.46, 5.96, 6.46 Xp. When a new subband is pulled below the Fermi 
level, the parallel Fermi sea built upon the newly occupied subband acquires more 
electrons, thereby increasing the screening and decreasing the energy-loss function. 
However, this effect is eventually overcome by the fact that all the subbands for the z 
motion get deeper with increasing film thickness. When a is increased by ~ Xp/2, a 
new subband begins to be filled and a new oscillation begins. The energy-loss function 
of a semi-infinite medium can then be extrapolated with the use of the following 



relation 20 



= lmW[a-}+lmW[a n } + lmWK} 
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where a~ = a n — A^/4 and a+ = a n + A^/4. 

The impact of short-range xc effects on the energy-loss function is investigated by 
employing the many-body kernel of Eq. (^|). These effects provoke a reduction in 
the screening of the electron-electron interaction, thereby increasing the energy loss 
for all projectile trajectories. This is observed in Fig. 2, where the RPA and ALDA 
energy-loss function is shown as a function of the z coordinate [z > in the vacuum, 
at the right-hand side of the slab] for two different values of a corresponding to a local 
minimum (a = an) and one of the two local maxima about the minimum (a = aJQj, 
and the same values of q\\ and u as in Fig. 1. The energy-loss function of the semi- 
infinite medium is also shown, by a solid line. We find that size effects become more 
significant as the particle moves away from the surface into the vacuum, which is due 
to the fact that as z increases the external field couples mainly with plasmon modes. 
Since these are modes characterized by their long wavelength, their behaviour depends 
strongly on the overall size of the system. For particles moving close to the surface the 
interaction has a short-range character, as the excitation of electron-hole pairs plays 
an increasing role. 

Fig. H] shows the RPA and ALDA probability P(uj), as obtained from Eq. (f|) 
versus the slab thickness a, for an external particle with v = 0.5 vq [vq is the Bohr 
velocity, v = 2.19 x 10 6 ms _1 ] and z = 2a to exchange energy to = 0.2 uj p with an Al 
slab. This figure exhibits an oscillatory behaviour similar to that shown in Fig. 1 for 
the energy-loss function. In Fig. |5] we have represented the RPA and ALDA energy- 
loss probability P(w), as a function of the energy transfer u, with the same values of 
v, z and a as in Fig. 3. As in the case of the energy- loss function, ALDA probabilities 
are well over those obtained in the RPA, due to the presence of short-range xc effects. 

The dependence of the RPA and ALDA stopping power on the thickness a of the 
slab, as obtained from Eq. (|IJ), is exhibited in Fig. ^. We have considered a charged 
particle moving with velocity v = 2vq along a definite trajectory at a fixed distance 
z = 2 ao from the right edge of an Al film. For these values of z and v the ALDA 
stopping power is significantly larger than that obtained in the RPA. Fig. shows the 
RPA and ALDA stopping power, as a function of the velocity, for a projectile moving 
at a fixed distance z = 2 from the surface into the vacuum. As the velocity increases 
the energy-loss spectrum of charged particles moving outside the solid is dominated 
by long-wavelength excitations and short-range xc effects, not included in the RPA, 
tend to become less important. 
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4 Conclusions 



We have investigated quantum-size effects on the energy-loss spectra of charged par- 
ticles moving parallel with metallic slabs, in the framework of linear-response theory. 

We have found that the quantization of the energy levels normal to the surface 
yields a neat oscillatory behaviour of the screened interaction, the energy-loss proba- 
bility and the stopping power, as a function of the thickness of the slab. The amplitude 
of these oscillations is found to decay approximately linearly with the slab thickness 
and their period is found to equal ~ We have also presented self-consistent 

calculations for a semi-infinite medium, which have been obtained from finte-slab cal- 
culations with the use of the extrapolation formula given in Ref. 20. 

Our results indicate that size effects become more significant as the particle moves 
away from the surface into the vacuum, where long-wavelength excitations play an 
important role. As for the impact of short-range xc effects on the various magnitudes 
that characterize the interaction of charged particles with metallic films, we have found 
that they increase the energy-loss probability, due to the reduction that these effects 
provoke in the screening of electron-electron interactions. 

A more detailed presentation of our self-consistent calculations of the energy-loss 
spectra of charged particles moving near a semi-infinite electron system will be pub- 
lished elsewhere|]15||. 
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Figure 1: The RPA and ALDA energy- loss function ImW^z, z; q\\,u), as a function of 
the slab thickness a, for q\\ = 0Aq F , to = 0.2 uj p , and z = 2a . Dashed lines represent 
the infinite-width limit, as obtained from Eq. @. The damping parameter in the 
evaluation of the density-response function of non-interacting Kohn-Sham electrons 
has been taken to be r\ = u p /10. 



Figure 2: The RPA and ALDA energy-loss function ImW^z, z; q\\,u), as a function of 
the z coordinate, for two different values of the slab width and the same values of q», 
uj and 7] as in Figure 1 . Thick solid lines represent the infinite- width limit of Eq. (^) . 
Dashed and thin-solid lines represent the energy-loss function for a = a u and a = a^, 
respectively. 



Figure 3: The RPA and ALDA energy-loss probability P{uj), as a function of the slab 
thickness a, for u = 0.2 u p , z = 2 a , v = 0.5 v , and r\ = u p /10. Dashed lines represent 
the infinite-width limit, which is obtained as in Eq. (||) with P instead of Imjy. 



Figure 4: The RPA and ALDA energy-loss probability P (u), as a function of the 
energy transfer oo, for two different values of the slab width and the same values of z, v 
and 7] as in Figure [5[ Thick solid lines represent the infinite- width limit. Dashed and 
thin-solid lines represent the energy-loss function for a = a n and a = a u , respectively. 



Figure 5: The RPA and ALDA stopping power, —(dE/dx), as a function of the slab 
thickness a, for z = 2a , v = 0.5 v , and r\ = u p /10. Dashed lines represent the 
infinite-width limit, which is obtained as in Eq. (|3|) with —(dE/dx) instead of ImVT. 



Figure 6: The RPA and ALDA stopping power, —(dE/dx), as a function of the velocity 
v , for two different values of the slab width and the same values of z and r\ as in 
Figure ||. Thick solid lines represent the infinite- width limit. Dashed and thin-solid 
lines represent the energy- loss function for a = an and a = a n , respectively. 
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